AC-field-induced Polarization for Uncharged Colloids in Salt Solution: 
A Dissipative Particle Dynamics Simulation 
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We study the response of a spherical colloid under alternating electric fields (AC-fields) by meso- 
scopic simulation method, accounting in full for hydrodynamic and electrostatic interactions. We 
focus on a special case of uncharged colloids. The main polarization mechanism is the "volume- 
polarization" , where ionic currents are deflected by the core of the uncharged colloid. Speciflcally, 
we compute the polarizability of a single colloid and systematically investigate the effect of AC 
field frequency and salt concentration. The simulation results are compared with predictions from 
classical Maxwell- Wagner theory and electrokinetic theory. 



I. INTRODUCTION 

Colloids suspended in an aqueous solution respond to 
external fields on relatively short time scales and in an 
often fully reversible way. This makes the application 
of electric fields one of the most attractive methods to 
manipulate single colloid or colloidal suspensions. Elec- 
trophoresis, the translation of individual colloid under a 
static external field is widely used to measure the 
surface charge density of colloidal particles. An alternat- 
ing electric field (AC-field) can also be used to probe the 
time-dependent properties of the colloidal suspension and 
provide substantially more information about the colloids 
than the static field. 

Let us consider a simple situation where a spherical 
colloid is immersed in either a dielectric or conducting 
medium such as a salt solution. The particle is often neg- 
atively charged, either by ionization or dissociation of a 
surface group, or by preferential adsorption of ions from 
the solution. The surface charges result in the formation 
of an electrical double layer (EDL) around the particle. 
The inner layer is called Stern layer and consists of ions 
which are absorbed on the particle surface. The outer 
layer is called diffuse layer and consists of counterions 
accumulated around the particle surface due to the elec- 
trostatic attraction between opposite charges. The ions 
in the diffuse layer are mobile and experience thermal 
motion. In Debye-Hiickel approximation, the thickness 
of the diffuse layer is characterized by the Debye length 
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where the sum runs over different ion types, Zi and piipo) 
are the valence and bulk concentration for i-th ion, re- 
spectively, and tm is the medium permittivity. 

When an electric field of the form E ~ Eoexp(ia;i) 
is applied to the suspension, both the colloidal particle 



and the electric double layer are polarized. The colloid 
acquires a dipole moment of the form p exp(ia;i), and the 
amplitude of the dipole moment can be written as 



p = a(a;)Eo, 



(2) 
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where a{ui) is the polarizability of the particle. The po- 
larizability is used to characterize the colloid response to 
the external fields, and it is a complex function of the fre- 
quency and the amplitude of the external field. For weak 
fields, the dipole moment is proportional to the magni- 
tude of the external field; thus the polarizability does not 
depend on the field strength. Various effects contribute 
to the dipole moment, which we will discuss in detail in 
the following. 

When no external fields are present, the diffuse layer 
of a charged colloid is spherically symmetric. Without 
loss of generality, we assume that the colloid is nega- 
tively charged (and hence the counterions are positively 
charged cations). The center of the counterion cloud and 
the colloid center coincide, resulting in a zero dipole mo- 
ment. When an external electric field is applied, the neg- 
atively charged colloid moves in the direction opposite to 
the field, while the positive counterions move in the di- 
rection of the field. The counterions move in and out of 
the electrical double layer, but when averaged over time, 
the counterion cloud has a dynamically stable ellipsoidal 
shape. The center of mass of the counterions is displaced 
with respect to the colloid center in the direction of the 
field, which results in a net dipole moment that points 
in the same direction. This type of polarization is re- 
ferred to as field-polarization. This mechanism can also 
be viewed as ion fluxes being deflected by the electric 
field produced by the charged colloid, as in ref. 

For field-polarization to take effect, the charged com- 
ponents should be mobile. The dipole moment is pro- 
duced due to the fact that oppositely charged compo- 
nents move in different directions in an external field. 
There exists another mechanism of polarization which 
can be traced down to the presence of the colloidal parti- 
cle as an obstacle that microions cannot penetrate. This 
mechanism has been described in the literature under 
different names, and for charged colloids, it is often re- 
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ferred to as concentration polarization 0, Fig. [Tfa) 
shows a schematic description of the mechanism. On 
the right-hand side of the colloid, a flux of negatively 
charged ions moves in the direction opposite to the field 
and hits the impenetrable core of the colloid. Combined 
with the abundant positive ions near the surface, the ac- 
cumulation of negative ions leads to an increase in the 
neutral electrolyte concentration. The slightly enhanced 
electrolytes produce a locally reduced Debye length, such 
that the counterions are compressed close to the surface. 
A similar process on the left-hand side leads to a decrease 
in electrolyte concentration and expansion of the coun- 
terion cloud. The net effect is a shift to the left of the 
center of the counterion cloud, and the resulting dipole 
moment points anti-parallel to the field. 



(a) (b) 




FIG. 1: (Colour on-line) The polarization mechanism 
due to the obstacle effect of the colloid core. Different 
interpretations exist for charged and uncharged colloids. 

(a) "concentration polarization" for a negatively 
charged colloid. The positive counterion cloud shifts to 
the left and the effective dipole moment points in the 
direction opposite to the external field. The bottom 
shows the electrolyte concentration profile ps along a 

line passing through the colloid center, (b) 
"volume-polarization" for an uncharged colloid. The 
bottom shows the charge density profile 
Pc = ^ P-); where p+ and p_ are the densities for 
cations and anions, respectively. 

The polarization mechanism due to the colloid obsta- 
cle is not restricted to charged particles. For uncharged 
particles, a dipole moment can also be induced by the 
movement of the salt ions. This is demonstrated in Fig. 
[IJb) . The positive ions move in the direction of the elec- 
tric field, and tend to accumulate on the left-hand side 
of the colloid due to the impenetrable core of the colloid. 
This produces a slight increase of the cation concentra- 
tion on the left-hand side of the colloid. A similar effect 
occurs on the right-hand side of the colloid where neg- 
ative ions accumulate, resulting in an increase of anion 
concentration. The combined effect again leads to an 
effective dipole moment which points in the direction op- 
posite to the external field. This is referred to as volume- 
polarization 0. In this variation of the charge 
density is produced near the colloid surface, whereas in 
the charged case, it is a variation of electrolyte density. 



The relative importance of these two types of polar- 
ization mechanism depends on the bare surface charges 
of the colloid, the bulk salt concentration, and the fre- 
quency and amplitude of the external field. It is easy to 
see that the field-polarization is only present when the 
colloids are charged. For charged particle, the situation 
is more complicated than for uncharged ones due to the 
presence of the electric double layer. Other effects also 
contribute to the polarization for charged particles, such 
as the lateral movement of microions in the Stern layer 
[1-01, and flow-induced effect due to the electrophoretic 
motion of the colloid d, . The situation is comparably 
simple for uncharged colloids, as the volume-polarization 
is the only mechanism for polarization. 

In this work, we study the case of uncharged colloids 
and focus on the effect of volume-polarization. This work 
is motivated by recent theoretical results due to Dhont 
and Kang 0], who gave accurate predictions for various 
response functions of uncharged colloids in salt solutions. 
We use Dissipative Particle Dynamics (DPD) simulations 
p^| - [T2| to study the dielectric response of a spherical col- 
loid under AC fields, including in full the hydrodynamic 
and electrostatic interactions. In a recent publication 
p^ . we have introduced a DPD model for colloids with 
no-slip surface, and presented preliminary results on the 
frequency-dependent dielectric response of nanoscale col- 
loids. The results were in surprisingly good agreement 
with the prediction of the Maxwell- Wagner theory, even 
though the latter was originally devised for micro-sized 
colloids. In the present work, we investigate the reason 
for this surprising agreement, and analyze in details the 
different physical factors that contribute to the polariz- 
ability of the colloid. To this end, we systematically vary 
the salt concentration, and we examine the contribution 
of the microions in different ion shells to the dielectric 
response. We compare our simulation results with a re- 
cent electrokinetic theory by Dhont and Kang Q and 
find very good agreement. The remainder of this article 
is organized as follows: in section |lT] we summarize the 
main theoretical results for uncharged colloids, which we 
will compare with our simulation results. In section Hill 
we give a brief introduction of the simulation model and 
describe important parameters of the system. We present 
the simulation results on the polarizability and effects of 
external AC fields in section IIVI Finally, section |V] con- 
cludes with a brief summary. 



II. THEORIES FOR UNCHARGED COLLOIDS 

In this section, we briefly review two types of theories 
for uncharged colloids: the Maxwell- Wagner theory and 
the theory by Dhont and Kang 0], which is based on 
the standard electrokinetic approach. We present only 
the important results in the main text and refer to the 
Appendices for detailed derivations. 

We consider a spherical particle of radius R, made of 
a material with permittivity Sp and conductivity Kp^ im- 
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mersed in a medium with permittivity and conductiv- 
ity In the simulation, we wih simphfy the picture by 
considering only non-conducting colloids (Kp ~ 0) and 
assuming that the particle and its surrounding medium 
have the same permittivity {cp = = e)- 

The classical Maxwell- Wagner theory [Tsl can be 
applied when there is a material property mismatch for 
the two components forming an interface. The condi- 
tion of continuity of the normal components of the dis- 
placement vector and the current density vector at the 
interface implies 



(3) 



If condition (jS]) is not satisfied, free ionic charges accu- 
mulate at the interface. The polarizability for a non- 
conducting colloid in a permittivity-matched medium 
(see Appendix |A| can be written as 



a{uj) = -~2neR 



with 



3 1 - ijuj/uJn-iw) 
1 + (w/cJmw)^ 



2K„ 
3e 



(4) 
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The Maxwell- Wagner time = l/co^yj corresponds 
to the time scale for the free charge formation at the 
interface. 

In the Maxwell- Wagner theory, the particle and its 
surrounding medium are treated as homogeneous sub- 
stances, and the induced charges only appear on the 
particle-fluid interface. The theory takes into account 
only the bulk properties of the components and ne- 
glects the effect of spacial distribution of the polarization 
charges. This simplification is valid for large particles, 
where the length scale for the variation is much less than 
the colloid size. For nanometer-size particles, the distri- 
bution of the microions near the interface may become 
important. 

The theory taking into consideration of the effect of 
space charge variation is based on the well-known, stan- 
dard electrokinetic equations Specifically, the theory 
for uncharged colloids describes the interplay of ion con- 
centration profiles and the electric potential. Two differ- 
ential equations are considered: the first one is the Pois- 
son equation, which connects the charge density to the 
Laplacian of the potential; the second one is the Nernst- 
Planck equation, which relates the time derivative of the 
ion concentration to the potential and flow velocity. A 
third equation, the Navier-Stokes equation, connects the 
flow velocity with the electrostatic force density. For the 
uncharged case, the convective contribution from the flow 
is a second-order effect and the Navier-Stokes equation 
can be omitted. The two equations (Poisson and Nernst- 
Planck) can then be solved with suitable boundary con- 
dition at the particle-fluid interface. 



In the literature, the majority of the studies have dealt 
with the general case of charged particles. For uncharged 
particles, Bonincontro et al. |16| obtained solutions for 
the electrokinetic equations and studied the effect of 
volume-polarization for a sphere in an 1-1 electrolyte 
solution. However, their result for the dipole moment 
differs from the Maxwell- Wa gner theory in the limit of 
kR — > oo. Later Garcia et al. 17[ derived a formula with 
the correct asymptotic behavior and showed that the ef- 
fect of space distribution of charges leads to a transition 
frequency higher than the Maxwell- Wagner results. Re- 
cently, Dhont and Kang have calculated the polarizability 
for an uncharged spherical colloid [2|: 



a{ui) 



1 + sR+ \ {sR)' 



s2 2 + 2si? + (si?)2 - i(Ki?)2 
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(see Appendix [B] for notations) . They also derived ex- 
plicit expressions of the response functions for the po- 
larization charge density and potential. Their results in- 
dicate that the characteristic length for the polarization 
charge variation is of the order of Id ■ 



III. SIMULATION MODEL 



In this section, we briefly review our simulation model 
for a colloidal particle in a salt solution and describe some 
important physical quantities. A more detailed descrip- 
tion can be found in ref. [ij . 

Our simulation system has three components: the sol- 
vents, the microions, and the colloidal particle. The sol- 
vent is modeled as a fluid of DPD beads, where DPD is 
used as a canonical thermostat (i.e. including the dis- 
sipative and stochastic part) without conservative forces 
[18| . In the following, physical quantities will be reported 
in a model unit system of a (length), m (mass), e (en- 
ergy), e(charge) and derived time unit r = a^Jmje. In 
these units, the temperature of the system is fc^T = 1.0 e, 
the density of the fluid is 3.0 cr^'^, and each solvent bead 
has a mass m. The DPD friction coefficient is set to 
IDPD = 5.0 sjmeja and the cutoff radius is Tc = 1.0a. 
The shear viscosity is measured using the method de- 
scribed in ref. jToj by analyzing a plane Poiscuille flow in 
a microchannel. For our system, we obtained a shear vis- 
cosity ris = 1-23 ± 0.01 m/r. which is in good agreement 
with ref. 0. 

Salt microions are introduced in the system as pairs 
of positively and negatively charged beads. We only con- 
sider the monovalent case where salt ions carry a single el- 
ementary charge ±e. The interactions between microions 
have two parts: one is the electrostatic interactions with 
the Bjerrum length Ib = e'^ / (ATremkBT) of the fluid set 
to 1.0(7. The other interaction is a short-range repulsive 
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Weeks-Chandler- Andersen interaction 12111 



VwcAir) 




i-^f + i for r<rc 
otherwise 

(7) 

whieh prevents the collapse of charged system. The cutoff 
radius is set at the potential minimum rc = vq + \f2a. 
The microions have a size of 1.0 cr (rg = 0). 

The colloidal particle is represented by a large sphere 
with interacting sites on its surface [23 - l27| . The short- 
range repulsive interaction between the colloid and fluid/ 
microion particles is modeled through a similar WCA 
potential ([7]). The radius of the colloid is i? = tq 4- 
a = 3.0 fT. To implement no-slip boundary conditions 
at the surface, a set of N interaction sites is distributed 
evenly on the surface. The interaction between these 
surface sites and the solvent beads is modeled using DPD, 
with the same DPD friction constant as that between 
fluid beads, but twice the cutoff range to ensure that a 
sufficient number of solvent particle can interact with the 
colloid. The mass of the colloidal particle is M = 100 m, 
and the moment of inertia is / = 360 mcr^, corresponding 
to a sphere with constant volume density. Fig. [5] shows 
a representative snapshot of a single colloidal particle in 
a salt solution. 




FIG. 2: (Colour on-line) Snapshot of an uncharged 
colloidal particle in a salt solution. The surface sites are 
represented by the blue beads. The red and green beads 
are cations and anions, respectively. Solvent beads are 
not shown here. 



We carried out simulations using the open source pack- 
age ESPResSo [1^. Modifications have been made to 
incorporate an external time-dependent electric field. A 
cubic simulation box of i = 30 cr with periodic boundary 
conditions in all three directions was used for all sim- 
ulations. Electrostatic interactions are calculated using 
Particle-Particle-Particle Mesh (P3M) method [H-iH. 
A time step of At = 0.01 t is used for the integration. 

As a benchmark test of our colloid model, we have per- 



formed simulations of an uncharged colloid in a salt-free 
solution. The electric field is turned off. We measured 
the mean-square displacement ((r(i) — r(0))^) of the col- 
loid as a function of time, shown in Fig. [31 Ballistic 
motion is observed for small times, as the mean-square 
displacement scales ~ t^. At late time, the colloidal par- 
ticle experienced many random collisions with the solvent 
beads, leading to a linear dependency on time. 
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FIG. 3: The mean-square displacement of an uncharged 
colloidal particle with radius i? = 3.0 ct in a salt-free 
solution. The linear regression is done for t > 30 r. 
Inset: normalized velocity autocorrelation function for 
the colloid, from ref. llSll. 



The diffusion constant of the colloid can be determined 
by a linear regression at late times. 



D = lim 



6t 



(8) 



The result after averaging five independent runs is D = 
0.014 ± 0.002 (T^/r. Alternatively, the diffusion constant 
can also be determined from the velocity autocorrelation 
function, using the Green-Kubo relation 



D 



1 



di(v(t).v(0)). 



(9) 



The velocity autocorrelation function has already been 
determined in our previous simulations [l3| and is shown 
in the inset of Fig. [S] The integration Eq. (j9l) gives 
the diffusion constant D = 0.013 ± 0.002 a^/r, in good 
agreement with the result from the mean-square displace- 
ment. For comparison, the diffusion constant of a Stokes 
sphere of radius R = 3.0(7 is D = ksT / {QTirjsR) = 
0.0144 ± 0.0002 crV'^: which indicates that the no-slip 
boundary condition is achieved in our colloid model. 

One important quantity required for comparison with 
the theory is the conductivity of the solution Km, which 
is related to the diffusion constant Di for microions 



2e^PsDi 



(10) 
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The formula is applied to 1-1 electrolyte solution and 
assumes that cations and anions have the same diffu- 
sion constant. The diffusion constant Dj can be de- 
termined by measuring the mean-square displacement 
for microions. We performed simulations with different 
salt concentrations, varying from ps = 0.003125 u^'^ to 
0.2 tj~^. We compared the simulation results with the 
empirical Kohlrausch law [s^ , which states that the mo- 
lar conductivity, or in our case Km/ Ps, depends linearly 
on the square root of the salt concentration ^fpl^ 



Kr, 

Ps 



knT 



= A-B^ 



'Ps 



(11) 



where A and B are fitting parameters. Fig. 0] shows the 
simulation results and a fit to Kohlrausch law. Note that 
the parameter A is equal to the limiting molar conduc- 
tivity [K„i/ps]p,-,o- 
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FIG. 4: The molar conductivity Km/ Ps, or the scaled 

diffusion constant of microions {ie^ /kBT)Di, as a 
function of salt concentration ps . The curve is a fit to 
Kohlrausch law with fitting parameters A — 1.50 and 
B = 2.15. 

In Table U we list the Debye length, the ratio R/Id 
between the colloid radius and the Debye length, the mi- 
croion diffusion constant, and the medium conductivity 
for four different salt concentrations. The simulations 
below are performed at these four salt concentrations. 

Ps [cT-^] Id [ct] kR Pi [^Vt] [eV(a£T)] 

0.003125 3.57 0.84 0.68 ± 0.03 0.0042 ± 0.0002 

0.0125 1.78 1.68 0.63 ± 0.04 0.016 ± 0.001 

0.05 0.89 3.36 0.50 ± 0.03 0.050 ± 0.003 



0.2 



0.45 6.73 0.27 ± 0.01 0.108 ± 0.004 



TABLE I: The Debye length Id, the ratio between the 
colloid radius and the Debye length kR, the microion 
diffusion constant Z?/, and the medium conductivity 
Km for different salt concentrations. 

To map the simulation units to real physical numbers. 



we use an aqueous solution of KCl as a reference sys- 
tem. The energy unit e is fc^T = 4.1 x 10~^^ J at 
room temperature T = 298 K. We set the Bjerrum 
length Ib = 1.0 cr, thus the length unit a corresponds 
to the Bjerrum length of water at room temperature, 
Ib = 0.71 nm. A salt concentration 0.2 cr"'^ in simula- 
tion translates to an experimental value of 912mol/L. 
We further equate the diffusion constant of microions 
at zero concentration {D'j = 0.75 (t^/t) to that of an 
electrolyte solution [s^. The diffusion constant for K+ 
and cr differ slightly (Dk+ = 1.96 x lO'^m^s-i and 
Dci- = 2.03 X lO-^m^s-i) [Mill, so we use the value 
2.0 X 10~^m^s~^ to obtain that one simulation time unit 
corresponds to real time 1.91 x 10^^" s. Thus the fre- 
quency / = 0.1 corresponds to an experimental fre- 
quency 523 MHz. 



IV. RESULTS AND DISCUSSION 

In this section, we report simulations for a single un- 
charged colloid of radius i? = 3.0 cr in solutions with dif- 
ferent salt concentrations. We systematically investigate 
the effect of varying the frequency of the external AC 
field, from / = 10"^ to 2.0 r~^. A separate simulation 
with a constant electric field is also performed to pro- 
vide the reference at the low- frequency limit. The am- 
plitude of the AC field is chosen in the linear region, 
Eq = 0.5 e/(cre). The dipole moment is calculated by 
counting the microions around the colloid. 



(12) 



where is the position of microion with charge qi. Since 
the total charge is zero, the dipole moment does not de- 
pend on the reference point and we use the colloid center 
Tc. We also need to choose a cutoff for the measure- 
ment and in this work we count the contribution from 
microions which have a distance less than 3 Id from the 
colloid surface. We have verified that the results do not 
change significantly when a larger cutoff is used. The 
time series of the dipole moment is collected during the 
simulation over 2 x 10^ time steps. The complex polariz- 
ability a{uj) is obtained by applying a Fourier transform 
to the time series of the dipole moment. 

Fig. [S] shows the simulation results for four differ- 
ent salt concentrations (see Table |l|, and correspond- 
ing Debye lengths ranging from 0.84 cr to 6.73 cr. Thus 
the Debye lengths are comparable to the colloidal radius 
{R — 3.0 cr). The real part of the polarizability Re{ci!} 
shows the in-phase component of the dipole moment with 
respect to the external AC field, while the imaginary part 
Im{a} gives the out-of-phase contribution. The polariz- 
ability in the low-frequency limit has a negative real part, 
indicating that the induced dipole is anti-parallel to the 
applied field, in agreement with the theoretical predic- 
tion. In the opposite limit of high frequency, the colloid 
and microions can no longer respond to the field, thus 
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FIG. 5: Real and imaginary components of the complex polarizability a{Lj) of an uncharged particle as a function of 
the frequency of AC field. The field strength is set in the linear region Eq = 0.5 e/ (ere). The salt concentrations are 
(a) ps = 0.003125(7^3^ (b) ^ 0.0125 cr-^, (c) p, = 0.05 cr'^, (d) p, = 0.2 a'^. The symbols are simulation results. 
The solid lines give the prediction from the Maxwell- Wagner theory with an effective radius i?cff = 2.92 cr. The 

dashed lines show the results of Dhont and Kang Q . 



both Re{a} and Im{a} converge to zero. 

The solid lines in Fig. [5] are predictions from Maxwell- 
Wagner theory. The limiting behavior of the polarizabil- 
ity at frequency / — >■ was obtained by an indepen- 
dent run in a DC field and used to calibrate the effective 
colloidal radius. Except for the case of the lowest salt 
concentration ps = 0.003125(7"'^ (Fig- EJa)): the low- 
frequency behavior for different salt concentration gives 
an effective radius ReS = 2.92 ± 0.07(7, which is reason- 
ably close to the "physical" radius R = 3.0(7. Having 
set the effective colloidal radius, we can calculate the 
Maxwell- Wagner prediction without further fit parame- 
ters. The resulting curves are in good agreement with 
the simulation data. Most notably, the theory predicts a 
crossover between two regions which is recovered in the 
simulation at roughly the predicted crossover frequency. 
The main discrepancy occurs at low salt concentration 
Ps = 0.003125 (7"'^ (Fig. EJa)), where the low-frequency 
limit of the Re{a} is higher than what the theory pre- 
dicts (with the assumption R^g = 2.92 a). This is prob- 
ably due to the finite size of the simulation box, as the 
Dcbye length for low salt concentration is quite large, 
and the resulting cutoff for calculating the dipole mo- 
ment is about half the size of the simulation box. In 



this case, the interaction between periodic images of the 
colloid may become important. We have qualitatively 
investigated the effect of finite size at higher salt concen- 
tration (smaller lo) by carrying out test runs of smaller 
systems. Finite size effects were found to affect the com- 
plex polarizability quantitatively, i.e., the value of Re(a) 
at low frequency and the peak value of Im(a). For linear 
system sizes L > 61 d, however, the transition frequency 
stays about the same. 

In the high frequency regime, the inertia effect of mi- 
croions and solvent also plays a role. When microions 
with finite mass are immersed in a viscous fluid, their re- 
sponse to the external field depends on the frequency. 
The characteristic time scale is tj = mj/Xi, where 
A/ = ksT / Di is the friction constant of microions. For 
solvents, there are two relevant time scales. One is the 
kinetic time ti, ~ o jv, defined as the time for momentum 
to diffuse over the characteristic length cr, and v is the 
kinetic viscosity of the solvent. Because the DPD sol- 
vent is compressible, the sound propagation will have an 
effect on the characteristic time scale of the sonic time 
^ca = ct/cs, where Cg is the speed of sound. We do not 
have a measurement value for the speed of sound, but 
from the velocity autocorrelation function (inset of Fig. 
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|3]) we estimate the time of acoustic momentum transport 
to be of order unity. In our simulation, ti ~ 3.7 r for 
Ps ~ 0.2a~^, and ^ 2.4 r. All these time scales are of 
order unity, therefore the system will be affected by in- 
ertia effects in the high frequency regime / > Ir^^. At 
low salt concentrations, the transition frequency is well 
below this value. At higher salt concentration, the tran- 
sition frequency increases and inertia effects come into 
play, leading to a slight reduction of the transition fre- 
quency. 

The dashed curves in Fig. [5] show the prediction of 
Dhont and Kang (also see Eq. ©). The difference 
between the two theoretic predictions is most prominent 
at low salt concentration, where the Dhont-Kang theory 
predicts a higher transition frequency than the Maxwell- 
Wagner theory. The polarization charges are distributed 
over a layer of width ^ Ijj (Debye length), and for small 
salt concentration, the Debye length becomes large. Thus 
the spatial variation of the polarization charges plays a 
more significant role when In ^ R. At the other limit 
of high salt concentration, the effect becomes less im- 
portant, and Eq. ^ can be shown to converge to the 
Maxwell- Wagner result, Eq. This can also be seen 
in Fig. Hfd) , where the two curves predicted by the differ- 
ent theories differ only slightly. It should be noted that 
the difference between the curves are smaller than the 
error of the simulation data at all salt concentrations. 
Thus the data on the polarizability cannot be used to 
assess the validity of the Dhont-Kang theory. 
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FIG. 6: The transition frequency /o as a function of the 
medium conductivity K^.- The transition frequency is 

calculated by fitting the imaginary part of the 
polarizability Im{a} by a Lorentzian. The solid and 
dashed lines are predictions from Maxwell- Wagner 
theory and Dhont-Kang results, respectively. 

Next we show the transition frequency /o of the dielec- 
tric response as a function of the medium conductivity 
Km in Fig. [5] The transition frequency is calculated 
by fitting the imaginary part of the polarizability Imja} 
by a Lorentzian. Also shown are the prediction from 
the Maxwell- Wagner theory (cf. Eq. (O) and numerical 
results based on the Dhont-Kang theory. When the ef- 



fect of polarization charge distribution is considered, the 
transition occurs at higher frequency than predicted by 
the Maxwell- Wagner theory. At high salt concentration 
(or large medium conductivity), the predictions of the 
two theories come close, but the transition frequency ob- 
tained from the simulations is lower than the prediction 
from both theories. The inertia effect mentioned above 
may contribute to this discrepancy due to the fact that 
/o approaches \t~^ for high salt concentrations. Fur- 
thermore, an analysis of the salt density distribution in- 
dicates that density oscillations appear near the surface 
for high salt concentration. In Fig. [3 the averaged ion 
densities are shown for two different salt concentrations. 
For low salt concentration {p^ = 0.0125 u"'^. Fig. El^a)), 
the ion density profile is monotonous and rises abruptly 
from near zero to the bulk value at the colloid surface 
(r < 3 <t). For high salt concentration (p^ — 0.2 a^^ , Fig. 
[71[b)), the density profile exhibits oscillations close to the 
surface. This layering effect of the microions near the 
colloidal surface at high salt densities is due to excluded 
volume effects between ions and may contribute to the 
reduced transition frequency in simulations. 
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FIG. 7: The cation and anion densities as a function of 
the distance to the colloid center. Densities for two 
different salt concentrations are shown here: (a) 
Ps = 0.0125 and (b) ps = 0.2 cr-^. The dashed lines 
show the bulk value. 

At low salt concentration (or small medium conductiv- 
ity), the Debye length is large compared to the colloidal 
radius, thus the effect of polarization charge is significant. 
The simulation results agree better with the Dhont-Kang 
result in this regime. The two theories differ in their as- 
sumptions how the polarization charges are distributed. 
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In the Maxweh- Wagner theory, the induced charges only 
appear at the interface, i.e. the colloidal surface, and a 
Laplace equation is assumed to hold in the solution re- 
gion. In the electrokinetic theory, a Poisson equation is 
applied in the solution region, and the Dhont-Kang re- 
sults predict a distribution of polarization charges with 
the characteristic length being the Debye length. Simu- 
lations can provide useful information about the charge 
density around the interface and verify the predictions 
from different theories. We shall demonstrate the ex- 
istence of spatially varying polarization charges by two 
methods. 

First, we can directly measure the local charge density 
of microions in the simulation. This is more easily im- 
plemented for a static electric field, as the charge distri- 
bution remains stationary and an average over long time 
can provide better signal-to-noise ratio. Fig. [5]shows the 
averaged charge density Pc{z) as a function of the dis- 
tance z from the center of the sphere along the direction 
of the external field. The curves are calculated using Eq. 
(|B12[) . The simulation demonstrates directly the spatial 
distribution of the polarization charges. There is an ex- 
cess of negatively charged ions on the front end of the 
colloid with respect to the external field, and the charge 
density decays to the bulk value over a distance charac- 
terized by the Debye length. For low salt concentration, 
the extent of the distribution is quite large, while for high 
salt concentration, it decays quickly to zero. The predic- 
tion from the electrokinetic theory and the simulation 
result agree quite well with each other. 

Secondly, we can divide the space surrounding the col- 
loid into spherical shells centered at the colloid and com- 
pute the contribution to the polarization from each shell. 
If the polarization charges are spatially distributed in- 
stead of being localized at the interface, the microions in 
the shell far away from the interface will contribute to 
the total polarizability. Using Eq. (jB8[) . we can compute 
this contribution from the electrokinetic theory. The re- 
sults for ps = 0.0125 cr~^ are shown in Fig. [5] for three 
different shells. The simulation results are in good agree- 
ment with the prediction, even the slight overshoot near 
the crossover region for the shell with large radius is re- 
produced. The contribution diminishes when we move 
away from the colloidal surface and reduces to noise level 
when the shell radius is more than three times of the 
Debye length. 



V. SUMMARY 

We have carried out mesoscopic molecular dynamics 
simulations of an uncharged colloidal particle in elec- 
trolyte solution under alternating electric fields, for dif- 
ferent salt concentrations. We have taken full account 
of the hydrodynamic interaction with thermal fiuctua- 
tions, using Dissipative Particle Dynamics, and the elec- 
trostatic interactions, using a Particle-Particle-Particle 
Mesh method. 



We obtained information about the dielectric response 
of the particle to a weak external AC field by comput- 
ing the complex polarizability of the colloid. We sys- 
tematically investigated the effect of the AC field fre- 
quency and the bulk salt concentration. The simula- 
tion results were compared with the predictions of the 
Maxwell- Wagner theory [l3, HH and from electrokinetic 
theory 0, Il6l [TtI . For the total polarizability, both the- 
ories show good agreement with the simulation, but the 
electrokinetic theory provides more information about 
the distribution of the polarization charges. 

There are some subtle differences between the elec- 
trokinetic model used in ref. and our simulation 
model. In ref. [H, the convective contribution in the 
Nernst-Planck equation is relatively small for uncharged 
colloids, thus the convective term is omitted. This 
is necessary to make analytical progress; it would be 
a formidable task to solve all equations including the 
Navier-Stokes equations. In the simulation model, the 
hydrodynamic interaction is included through the dissi- 
pative particle dynamics. The good agreement between 
the simulation and theory indicates that, for uncharged 
colloids, the hydrodynamics is not important. However, 
this is not true for charged colloids, where the contribu- 
tion from the electrophoretic motion of the colloid and 
the electro-osmotic flow alter the results significantly. 

Secondly, the salt microions are considered to be point- 
like particles in the electrokinetic equations, which are 
introduced through the ionic concentration. This treat- 
ment corresponds to a mean-field type approach and ne- 
glect the ion-ion correlations. The effect of ion-ion corre- 
lations is small for the system of low salt concentration, 
but becomes important at high salt concentration. In 
the simulation, we used a particle-based approach, and 
microions are modeled as spherical beads with finite size. 
As a consequence, we observed the layering effect of the 
microions near the colloid surface, which cannot be de- 
scribed by a theory that neglect ion-ion correlations. It 
will be interesting to compare the simulation results to 
the electrokinetic model which includes the ion size effect 

[Mill. 
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Appendix A: Maxwell- Wagner theory 

In this appendix, we give a short introduction to 
Maxwell- Wagner theory. More detailed information can 
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FIG. 8: The charge density pc = e(p+ — p ) as a function of the distance from the colloid center along the direction 
of the external field. The charge densities have been normalized using the bulk salt concentration ps . The results for 
four different salt concentrations are shown here (see Table |T|). The curves are predictions from Eq. (jB12l) . 
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FIG. 9: Contribution to the polarizability from microions located in a spherical shell at distances (a) r S [3,4], (b) 
[4,5], and (c) [5,6], respectively. The salt concentration is ps = 0.0125 cr""^ and the corresponding Debye length is 
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be found in refs. 



When loss is present, the dipole moment of a particle 
immersed in fluids exhibits a phase lag with respect to the 
external AC field Eq exp(iwt). A complex dipole moment 
can be written as 



(Al) 



where the Clausius-Mossotti factor K{u!) is a complex 
number containing both the magnitude and the phase 
information about the effective dipole moment. In the 



Maxwell- Wagner theory, it has the form 

f* — fl^ 

K(e* e* ) 



(A2) 



where e* and are the complex dielectric constants of 
the particle and the medium, respectively. They are de- 
fined as 



iuj 



e„, + —, (A3) 



where e (without the star) and K are the permittivity and 
conductivity, respectively. The real and imaginary parts 
of the complex Clausius-Mossotti factor can be expressed 
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Re{K*} 



Kp ~ Km 1 + CJTQTmw 
Kp + 2Krn 1 + CjV2„ ' 



(A4) 



Kp + 2Km 1 + wV2^ ' 

6r) ~ ^m. 



K — K ' 



(A6) 
(A7) 



For e™ = Ep = e and Ap — 0, the Clausius-Mossotti 
factor is reduced to 



Re{A:*} 



2 1+ (cj/Wniw)^ 



(A8) 



(A9) 



2 1 + (a;/a;n,w)2 ' 
with Wniw = 2A'm/(3e). The induced dipole moment is 

3 1 - i(a;/w,nw) 



(AlO) 



1 + (aj/a;,„w)2' 
This corresponds to the polarizabihty given in Eq. (j4|. 

Appendix B: Dhont-Kang theory 

Dhont and Kang Q derived a formula for the com- 
plex response function of the polarization charge density. 
Here for the purpose of comparison with simulations, we 
derive two equations (|B8p and (|B12p for our system from 
their results. The system consists of an uncharged non- 
conducting colloid immersed in an electrolyte solution 
with the same permittivity as the colloid core. The ex- 
ternal electric field has a form of E = Eq exp(iajt) and 
the charge density can then be written as 



Pc{r,t) = R{r,Lu)exp{iujt), 



(Bl) 



where R is the complex- valued response function for the 
polarization charge density. The response function takes 
the form 



i?(r,cj) = i?o(w)Eo • V 



exp(— sr) 



with 



en R-^ exp(si?) 



2 + 2si? + (si?)2 - i(Ki?)2' 

k{x + iy), 



y 

A 



1 

V2 
1 

UJ 

K^Dl' 



1 + (1 + A2)1/2' 



1/2 



-1+(1+A2)1/2 



1/2 



(B2) 

(B3) 
(B4) 
(B5) 

(B6) 
(B7) 



From the response function R, the contribution to the 
dipole moment from the shell ri < ?' < r2 can be readily 
calculated 



P(w|ri,r2) 



dr ri?(r, ut) 



. ^ ^ , ,exp(— sr) 



1 2 

1 + sr + - (sr) 



(B8) 



The total dipole moment is then obtained by taking ri 
R and r2 — > oo, 



p(w) = p{uj\R,(X)) 



-47reEoi?' 



k2 l + .si?+i(si?)2 
^ 2 + 2si? + (.si?)2 - 



(B9) 



This is Eq. (30) in ref. 0, and it can be shown to be in 
agreement with Eq. (19) in ref. [l7j . 



Eq. ((B9| can be rewritten using (|B4| - ((B7| and taking 
the limit sR — > oo. 



p(w) = -27reEoi?^ 



1 + |Ai' 



(BIO) 



The factor |A in the denominator can be expressed in 
terms of AT,,,, using Eq. ^j}, (UHl) and ([T|), 



3 . 3ll) Swe UJ 

-A = — ^ = = 

2 2k^_D/ 2A',„ ujmv 



(Bll) 



This is exactly the result from the Maxwell- Wagner the- 
ory (cf. Eq. (|JTO)) ). 



For a constant electric field, the imaginary part of the 
response function vanishes. The real part of R, which 
in this case (w = 0) is just the charge density pc, as a 
function of the distance from the center of the colloid 
along the direction of the external field is 



p,{z) = Re{R{z,Lu = 0)} 

fi^R-^ 1 + Kr 



2 + 2KR+l{nRY 



-K{r-R) 



(B12) 
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